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ABSTRACT 



The cosmic distance scale largely depends on distance determinations to Local Group galaxies. In this sense, the Andromeda Galaxy 
(M 3 1) is a key rung to better constrain the cosmic distance ladder. A project was started in 1999 to firmly establish a direct and accurate 
distance to M 3 1 using eclipsing binaries (EBs). After the determination of the first direct distance to M 3 1 from EBs, the second direct 
distance to an EB system is presented: M31V J00443610+4129194. Light and radial velocity curves were obtained and fitted to derive 
the masses and radii of the components. The acquired spectra were combined and disentangled to determine the temperature of the 
components. The analysis of the studied EB resulted in a distance determination to M 3 1 of (m - M)o = 24.30 ± 0. 1 1 mag. This result, 
when combined with the previous distance determination to M31, results in a distance modulus of (m - M)o = 24.36 ± 0.08 mag 
(744 ± 33 kpc), fully compatible with other distance determinations to M31. With an error of only 4%, the obtained value firmly 
establishes the distance to this important galaxy and represents the fulfillment of the main goal of our project. 
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1. Introduction 



Stars: fundamental parameters - Stars: distances - Cosmology: distance scale - Galaxies: 



Eclipsing binaries (EBs) have always been an important t ool for 
testing and determin i ng the physical prop erties of stars (Popper 
Il967t iGuinanl 11993b iTorres et ail 120091) . They are composed 
of two stars that, when orbiting each other, produce periodic 
eclipses. The great potential of EBs is that their orbital mo- 
tion, inferred from the radial velocity curves, and the shape 
of eclipses, obtained from the light curves, can be entirely ex- 
plai ned by the grav itation laws and the geometry of the system 
! (see lHilditchll2001l for details). 

The direct determination of the radii of the components 
of EB systems mad e that several authors (e.g., lLacvl Il977t 

iGimenez et all 1994 suggested the possibility of using EBs for The main interest for an accurate distance determination to 



is the surface brightness or, equivalently, the effective tempera- 
ture of the components. 

The potential of EBs to derive distances encouraged sev- 
eral projects to obtain dire ct distance determ inations, either 
within the Mil ky Way (e.g..lMunari et al]|2004l) . the Magellanic 
Clouds (e.g. jFitzpatrick et al.ll2003|, and references therein), or 
M31/M33 (DIRECT project. iKaluznv et al.lll998l) . In 1999, a 
new project was started to obtain a direct dista nce determina- 
tion to the Andromeda Galaxy (M 31) from EBs (Ri bas & Jordil 
120031; iRibasetal. .120041) , prov iding the fi rst direc t distance deter- 
mination in Rib as et all {2005, hereafter Pa per J) . 



deriving distances. The only additional requirement to determine 
the absolute luminosity of an EB system and, hence the distance, 



* Based on observations made with the Isaac Newton Telescope op- 
erated on the island of La Palma by the Isaac Newton Group in the 
Spanish Observatorio del Roque de los Muchachos of the Instituto de 
Astroffsica de Canarias. 

** Based on observations obtained at the Gemini Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under a cooperative agreement with the NSF on be- 
half of the Gemini partnership: the National Science Foundation (U- 
nited States), the Science and Technology Facilities Council (United 
Kingdom), the National Research Council (Canada), CONICYT (Chil- 
e), the Australian Research Council (Australia), Ministerio da Ciencia e 
Tecnologia (Brazil) and Ministerio de Ciencia, Tecnologfa e Innovation 
Productiva (Argentina) 



M 3 1 lies on the pote ntial of this galaxy to b e a first-class dis- 
tance calibrator (e.g.. IClementini et al.l l2Q0lh . The reasons for 
this are: (1) Contrary to the Magellanic Clouds, the distance to 
M31 is large enough so that its geometry does not introduce 
any systematics in the final distance determination; (2) typically 
with a moderate re ddening value (E(B - V) = 0.16 + 0.01, 
Mass ev et al.l ll995). it is close enough to enable the individual 
identification of stars suitable for distance determination (such 
as EBs or Cepheids); (3) an Sb I— II giant spiral galaxy (like 
M 31) provides an appropriate local counterpart for the galaxies 
comm only used for distance determination (e.g jFreedman et all 
2001); and (4) M31 can also provide an absolute calibration of 
the Tully-Fisher relationship, enabling the calibration of the fur- 
thest distance determination methods. Therefore, the character- 
istics of this spiral galaxy make it an important step of the cosmic 
distance scale. 
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As mentioned above, our project al ready pr ovided the first 
direct distance determination to M31 dPaper J) . However, the 
excellent dataset obtained (Sect. [2j allowed the determination of 
additional distances in order to further constrain the distance to 
M 3 1 . Therefore, in the present work, we are presenting the sec- 
ond direct distance determination to M 3 1 from an EB system 
(Sects. [3|-[6]l. This result (the last one until further spectroscopic 
observations can be secured) enables a critical discussion on the 
distance to M31 (Sect.[7]i. 

2. Observations 

Photometric time series (in B and V passbands) were acquired 
with the Wide Field Camera at the 2.5 m Isaac Newton Telescope 
at La Palma (Spain). The field of observation covers 34' x 34' 
at the North-Eastern part of M 3 1 . Around 260 images were ob- 
tained in each passband during five campaigns between 1999 and 
2003 . The Difference Image Analysis technique (DIA. IWozniakl 
2000) was used to per form the photometric data reduc tion (thor - 
oughly explained in Vilard ell et al.l [2006, hereafter iPaper III) , 
providing light curves for more than 3 964 variable stars, with 
437 being identified as EBs. In add ition, the V photo metric 
observations of the DIRECT project (Stan ek et al.l ll999) were 
also considered. Although the observations are somewhat nois- 
ier, they allowed us to extend the time baseline to 1996 and to 
ensure the consistency of our DIA photometry. 

The EB sample was inspected to d etect tho se targets more 
suitable for distance determination (see Pa per It) . A 5' x 5' re- 
gion containing five of the selected EBs was observed with the 
multi-object spectrograph (GMOS) at the Gemini-North tele- 
scope (program ID GN-2004B-Q-9). The instrumental setup was 
set to provide the highest possible resolution (R=3 744) using 
the slit width of 0'.'5 and the B1200_G5301 grism. The spectra 
cover the wavelength interval between 390 and 530 nm, with two 
gaps between 436.7-436.9 and 485.8-487.0 nm. A total of eight 
exposures were taken between September and November 2004 
with an integration time of 4100 seconds each and an additional 
one with an exposure time of 3240 seconds in February 2005 
(Table [TJ. The data reduction was performed with the Gemini 
IRAF package version 1.7, providing spectra with a S/N between 
6 and 39 for the five selected EBs. 

After careful analysis of the five EBs, two of the observed 
targets resulted to be optimum for obtaining a precise dis- 
tance determination to M31: M31V J00443799+4129236 and 
M31VJ 004436 10+4129194. The first target was already ana- 
lyzed in lPaper j and the analysis of the second target is presented 
here. The three remaining EBs are unsuitable for distance deter- 
mination. One of them is a single-line EB and, therefore, the 
absolute properties of the stars cannot be obtained. Another EB 
is too faint for an accurate analysis of its spectra. The third one 
is a triple-line system, where the important contribution of the 
third light prevents a precise distance determination. 



3. Radial velocities 

The determina tion of radial velocities ( RVs) was performed, as 



in the case of lPaperll with TOD COR (IZucker & Mazehl 



faper H, with lUU tUK (IZu cKer & Mazer 
and the ATLASSQ and TLUST\0 ( Lanz & Hubenvll2003 



1994) 



2007) 



synthetic models. In this case, the best pair of synthetic spec- 
tra was determined with t he following iterative appro ach. As a 
first step, the preliminary IWilson & Devinnevl dl97ll hereafter 



1 Available at: http : //kurucz . cf a . harvard . edu/ 

2 Available at: http://nova.astro.umd.edu/ 



IW&Dh fit (Sect.H was used to define an initial list with pairs of 
models having a temperature ratio, g ravity ratio and rotational 
velocities compatible with the IW&Di parameters. All the model 
pairs in the initial list were then used to determine RVs and a 
simple RV curve model was fitted to the obtained values. The 
free parameters in the fit were the semi-major axis (a), the sys- 
temic velocity (y) and the mass r atio (q) , while the period and 
reference time was fixed from the IW&Dl solution. From all the 
derived solutions, the models having a lower dispersion around 
the fitted RV curve were selected and all the neighboring models 
in the space of parameters were also attempted. The process was 
repeated until the local minimum was found, (i.e., none of the 
neighboring models has a scatter lower than the selected pair of 
models). 

Since the resulting RVs depend, to some extent, on fhe lW&Dl 
fit, the obtained values were us ed to fi nd a new solution with the 
IW&Dl (Sect. |H. The resulting IW&Dl solution was then used to 
determine new RVs. The process was repeated until the pair of 
synthetic models providing the best fit was the same in two iter- 
ations. The final solution shown in Table [T] contains all the RVs 
(corre cted to the heliocentric reference frame) that remained af- 
ter the W&D fit (with 3cr clipping). Rejected observations shown 
in Fig.Q]were obtained during eclipses (one of them corresponds 
to the spectrum of February 2005 with a shorter exposure time). 

The rejection of the observation at phase 0.885 cannot be ex- 
plained neither by the proximity to the nodes (since other obser- 
vations are obtained at larger phases) nor by a lower S/N of the 
observed spectrum. After some additional tests, it was observed 
that the rejected observation could be recovered with different 
pairs of templates at the cost of losing other observations close 
to the nodes and a larger dispersion of the fit. In all the different 
solutions attempted, the systemic velocities and semi-amplitudes 
were in perfect agreement with the solution presented here. 

4. Mass and radius determination 

The mass and radius determin ation wa s performed following the 
fitting procedure described in lPaper j . Both semi-detached con- 
figurations (with either the primary or the secondary filling the 
Roche lobe) yield fits with the same residuals. The final configu- 
ration was adopted after the temperature determination analysis 
(Sect. [5]). The spectra resulting from the disentangling indicate 
that the secondary is much fainter than the primary, while the 
assumption of the secondary filling its Roche lobe yields a flux 
ratio of Fy,s/Fv,p = 0.85. Therefore, a configuration where the 
primary (instead of the secondary) fills the Roche lobe was used 
in this case. 

In order to ensure the viability of the adopted scenario, and to 
rule out any other possible configurations, additional fits adopt- 
ing a detached configuration were performed. The resulting pa- 
rameters revealed that, depending on the initial value of the sur- 
face potential of the secondary component, either the primary or 
the secondary tended to fill their respective Roche lobes. 

With the adopted configuration, the final rms residuals are 
0.014 mag in B, 0.015 mag in V, and 0.047 mag in the DIRECT 
V light curve. The residuals of the RVs are 13 and 6 km s _1 
for the primary and secondary components, respectively. The 
light and RV curves, with their respective fits superimposed, 
are shown in Fig. Q] The resulting best-fitting elements listed 
in Table [2] reveal two components with masses and radii of 
M P = 21.7 ± 1.7 M , Rp = 9.2 + 0.2 R and M s = 15.4 + 1.2 
M Q ,R S =5.6 + O.4R . 

Light curves measured using DIA photometry require an ac- 
curate reference flux determination to avoid any bias in the scale 
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Date 


Time 


Exp. time 


S/N 


HJD 


Phase 


Primary 


Secondary 




[UT] 


LsJ 








[km s~'] 


[km s- 1 ] 


Nov 7, 2004 


10:14 


4100 


22 


2453 316.931 


0.272 


-370.1±8.2 


143.5+30.3 


Sep 15, 2004 


13:29 


4100 


11 


2453 264.066 


0.467 






Feb 12, 2005 


5:39 


3240 


8 


2453 413.734 


0.524 






Nov 12, 2004 


7:35 


4100 


15 


2453 321.820 


0.659 


14.3+13.9 


-414.1+25.5 


Nov 12, 2004 


8:49 


4100 


16 


2453 321.872 


0.684 


57.3+13.5 


-448.2+24.8 


Nov 10, 2004 


9:00 


4100 


24 


2453 319.879 


0.711 


30.6+9.9 


-455.5+17.4 


Sep 14, 2004 


8:53 


4100 


20 


2453 262.874 


0.885 






Sep 12, 2004 


9:45 


4100 


13 


2453 260.910 


0.927 


-69.7+18.1 


-309.2+19.4 


Oct 17, 2004 


6:31 


4100 


14 


2453 295.776 


0.946 


-58.1+24.2 


-271.7+22.9 



Table 2. Fundamental properti es of M 3 1 V J004436 1 0+4 1 29 194 
derived from the analysis with W&Dl 



System properties 


Ca* (mag) 


19.832+0.013 


VSm (mag) 


19.948+0.015 


P (days) 


2.048644+0.000007 


t mm (HJD) 


2452 908.694+0.004 


i (deg) 


69.9+1.6 


y (km s _I ) 


-164+5 


a (R e ) 


22.6+0.5 


q = Ms /Mp 


0.71+0.04 


T c ff,s /T c ft\p 


0.897+0.019 


Fb,s 1 Fgp 


0.33+0.03 


Fv,s 1 Fyp 


0.33+0.03 


Fd,s 1 F" DP 


0.33+0.05 


Component properties 


Primary Secondary 


k(Rq) 


9.2+0.2 5.6+0.4 


M(M Q ) 


21.7+1.7 15.4+1.2 


log g (cgs) 


3.85+0.02 4.12+0.05 


K b (km s -1 ) 


210+9 306+11 


v sync sin / (km s _1 ) 


213+6 131+8 



" Out of eclipse average: = [0.20 - 0.30, 0.70 - 0.80] 



Including non-Keplerian corrections 



of the light cu rves. Following the same procedure used for the 
EB in iPaper J, several light-curve fits were performed with a 
variable third-light contribution (It,). Again, the excellent agree- 
ment with DIRECT light curve and the convergence of the fits 
to I3 ~ 0, ensures that the flux zero-point is correct and rules out 
any possible relevant blends with unresolved companions. 



Table 3. Fundamental properties of M3 1 V J004436 1 0+4 1 29 1 94 
and resulting distance determination, derived from the model- 
ing of the optical spectrum of the primary component, using 
TLUSTY atmosphere models and the values in Table|2] 



System properties 


M v (mag) 


-4.90= 


=0.08 


E(B - V) (mag) 


0.18d 


=0.02 


A v (mag) 


0.55= 


=0.08 


(in - M) a (mag) 


24.30= 


=0.11 


Component properties 


Primary 


Secondary 


7W (K) 


33 600+600 


30 100+900 


logg (cgs) 


3.86+0.12 




v rot sin; (km s -1 ) 


189+12 




My (mag) 


-4.59+0.07 


-3.38+0.12 


(B - V)o (mag) 


-0.295+0.002 


-0.286+0.004 



ity. The resulting values are listed in Table [3] together with the 
values of absolute luminosity (M v ) and intrinsic color ((B-V)o). 
These values wer e computed by scaling the absolute magnitude 
in the models of Lejeune & Schaerer (2001) (for a given tem- 
perature and gravity) with the observ ed radii of the stars. When 
needed, the values resulting from the IW&Dl fit were used (e.g., 
to determine the temperature of the seco ndary). This procedure 
is equivalent to the one used in Paper I, where instead of per- 
forming the synthetic photometry to obtain the surface fluxes, 
they were obtained from the models. 

Finally, once the absolute magnitudes of the components are 
determined, the determination of the distance is straightforward. 
The resulting distance modulus to M31V J00443610+4129194 
is (m-M) = 24.30+0.1 1 mag or, equivalently, d = 124+37 kpc. 



5. Temperature and distance determination 

The last required ingredient to obtain a direct distance deter- 
mination to M31 is the temperature and line-of-sight absorp- 
tion det ermination. Again, following the same procedure as in 
Paper I, we determined the temperature of the components by 
modeling their disentangled spectra. The s pectral disentan gling 
was performed with the KOREL program (Hadrava 1995) with 
all the req uired parameters fixed to the values determined from 
the IW&Dl analysis (Sect.g). 

Considering the low S/N of the spectrum of the secondary, 
the fitting procedure was restricted to the primary compo- 
nent. The temperature (T e ff), surface gravity (logg) and rota- 
tional velocity (v lot sin 2) were fitted considering TLUSTY tem- 
plates. Spectra with solar metallicity were adopted, since they 
were found to provide good fit s to th e observed spectrum. 
In addition, the EB analyzed in Paper I, only 29 arcsec from 
M31V J004436 10+4 129 194, was found to have solar metallic- 



6. Comparison with stellar evolutionary models 

The final step towards the characterization of 
M31VJ004436 10+4 129 194 is the comparison of the de- 
rived physical properties with stellar evolutionary models. The 
comp arison has been ma i nly p erformed with the Geneva models 
of iLeieune & Schaerer! d2001l) . con sidering sola r metallicity. 
Other models have been considered dClaretll2004l) without any 
major variations on the derived conclusions. It is important to 
emphasize that the evolutionary tracks are built for isolated 
stars. In close binary systems, the evolution of both stars can 
be largely modified, with respect to their isolated evolutionary 
tracks, when one of the components fi lls the Roch e lobe and 
mass transfer takes place (see, e.g., IVanbeverenl 1 1993b . In 
addition, due to the effects of mass loss in massive stars (due to 
stellar wind), the current derived masses of the components will 
be different from the masses at the Zero-Age Main Sequence 
(ZAMS). Therefore, we initially fitted different evolutionary 
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Fig. 2. Comparison of the individual disentangled spectrum for the primary component of M31V J004436 10+4 129 194 (above) 
with the TLUSTY synthetic spectrum (middle). The bottom spectrum shows the noisy disentangled spectrum of the secondary 
component. 



tracks in the mass-radius diagram to derive the mass of each 
component at the ZAMS (AWams = 22.1 M Q , Ms zams = 15.4 
Mo). 

Once the evolutionary tracks are adopted, the Hertzsprung- 
Russell (H-R) diagram can be studied. The location of the com- 
ponents on the H-R diagram (Pig. [3j reveals that both agree with 
their predicted evolutionary tracks. The derived properties seem 
characteristi c of a detached system. However, as explained in 
Sect. H] the IW&Dl solutions supposing a detached configura- 
tion converge to solutions where either the primary or the sec- 
ondary fill their respective Roche lobes. In addition, the dis- 
entangled spectra seem to favor the case where the primary is 
filling the Roche lobe, since it is much brighter than the sec- 
ondary. All these observations can be explained supposing that 
M31VJ00443610+4129194 is a pre-mass transfer EB system, 
where the primary component is almost filling its Roche lobe. 

In a pre-mass transfer EB, both components are detached and 
basically should follow the evolution of single stars. This ex- 
plains why the observed properties agree with the evolutionary 
tracks, since each component has evolved without any major in- 
teraction with the other component. This could also explain why 
the most massive component can be assumed to be filling the 
Roche lobe. What is generally observed in semi-detached EBs 
is the Algol paradox, where the do nor is th e less massive com- 
ponent (as in the case of the EB in Paper I). The reason of this 
phenomenon dPaczvriskill 19711) is that, when the donor is more 
massive than the companion, the mass transfer is very rapid, tak- 
ing place on a thermal time scale (~ 10 5 yr). On the contrary, 
when the donor is the less massive companion, the mass transfer 
can last as long as the main-sequence life-time of the donor. The 
lack of any signature of intense mass transfer (such an O'Connell 
effect in the light curves), combined with the short timescale of 
the process, makes that a situation where the most massive com- 
ponent is filling the Roche lobe is unlikely. However, the time 
spent in a situation where the massive component can be almost 
filling the Roche lobe has a typical timescale of ~ 0.5 Myr, mak- 
ing the situation more likely to be observed. 



The close agreement of the observations with the theo- 
retical models, apart from clarifying the evolutionary stage 
of M31VJ004436 10+4 129 194, also allowed the fitting of an 
isochrone, resulting in an estimated coeval age of 4.2 + 0.4 Myr. 



7. The distance to M31 

The two analyzed EBs can be used to provide a distan ce deter- 
mination to the M 31 galaxy (Table|4j>. As mentioned in lPaper j, 
the derived distances to each EB can be associated to the center 
of M31, because the correction due to the location of the EB 
within the galaxy is negligible (~ 0.3%). Therefore, the two in- 
dependently derived distances, in mutual agreement within their 
one sigma error bars, show that EBs can be used to derive precise 
and accurate distances to M 3 1 . 

The two double-line EBs have been combined to derive a 
weighted mean distance to M 3 1 of 744 + 33 kpc or (m - M)o 
24 .}() + 0.08 mag. This resulting distance relies on the modeling 
of two different EBs and, therefore, can be considered to be: 

Direct. The distance determination from EBs does not rely on 
previous calibrations and, since it is derived using a one-step 
procedure, it can be considered direct. Therefore, a recalibration 
of the distance to other Local Group galaxies (such as LMC) 
or a variation on the zero-point of any standard candle (such as 
Cepheids) has no effect on the derived distance. Furthermore, 
any standard candle in M 3 1 can be calibrated using our derived 
distance. 



Precise. The uncertainty on the derived distance modulus of 
0.08 mag represents a distance determination with an error of 
only 4%, which is remarkable, given the faintness of the stud- 
ied targets. Although other distance determina tions to M 3 1 with 
smaller errors have been reported (see a list in lPaper IJ) . most of 
these results do not consider the effect of systematics. In fact, the 
derived uncertainty is equal to the standard deviation resulting 
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EB system 


V 


M v 


A v 


(m - M) 


Distance 




[mag] 


[mag] 


[mag] 


[mag] 


[kpc] 


M3 1 V J00443799+4 1 29236 


19.27+0.02 


-5.77+0.06 


0.60+0.10 


24.44+0.12 


772+44 


M3 1 V J004436 1 0+4 1 29 1 94 


19.948+0.015 


-4.90+0.08 


0.55+0.08 


24.30+0.11 


724+37 


Weighted mean 








24.36+0.08 


744+33 



fr om the c ombination of the non-direct distance determinations 
in lPaper III ((m - M)o = 24.39 + 0.08 mag) and, therefore, our de- 
rived distance is equally precise. The uncertainty in the derived 
distance could be improved with the analysis of addit ional EB s 
(already identified in the list of 24 EBs provided in iPaper III) , 
resulting in a distance determination to M3 1 with a relative un- 
certainty of 2-3% and free of most systematic errors. This result 
would represent the most accurate and reliable distance determi- 
nation to this important Local Group galaxy. 

Accurate. One of the most important points when determin- 
ing distances is the effect of possible systematics in the derived 
value. Contrary to other distance determinations, the uncertainty 
in our distance modulus includes most, if not all, the possible 
systematics. In particular, the possible sources of systematic er- 
rors, and the corresponding considerations, can be summarized 
in the following points: 

- Photome try. The photometry has been compared with other 
catalogs (IPaper Ij) and checked to be well below 0.03 mag 
for the magnitude of the selected EBs. 

- Assumed configuration in the modeling of the EBs. The con- 
figuration assumed is completely independent for each one 
of the two EBs used for distance determination and results 
in distances that agree within the uncertainties. In addition, 
M31VJ00443799+4129236 has clear evidences of being a 
semi-detached EB (O'Connell effect, no eccentricity, etc.) 
and M31V J004436 10+4 129 194 has been thoroughly tested 
for any other possible configuration, with none of them being 
capable to reproduce the observations. 

It has traditionally been argued that detached EBs are the 
only systems capable to provide accurate distances. The 
most common reason is that non-detached EBs are affected 
by the proximity of the components, introducing distortions 
and reflection effects. However, the proximity of the com- 
ponents can properly be t aken int o account by current mod- 
eling algorithms (such as IW&Dh . In addition, the fact that 
one of the components fills the Roche lobe decreases the 
number of free para meters and greatly impr oves the stabil- 
ity of the solution (IWvithe & Wilson! [2002). On the other 
hand, there are some other effects, missing in detached EB, 
that could introduce some systematics in the solutions (e.g., 
hot spots, circumstellar disks, etc.). In any case, these ef- 
fects can be observed from the acquired data and prop- 
erly modeled when they have a mild effect on photome- 
try. A clear example of these effects is the presence of the 
O'Connell effect in M31 V J00443799+4129236. The mod- 
eling performed, with the introduction of a hot spot, enabled 
the use of M31V J00443799+4129236 for distance deter- 
mination. Therefore, using semi-detached EBs reduces the 
number of free parameters, at the cost of having to account 
for some other effects that, when properly considered, should 
introduce no systematic error in the derived distances. 

- Radial velocities. Several aspects have been considered with 
respect to the radial velocity determinations. TODCOR is a 



well tested program that has been checked to introduce some 
syste matics only for spect ra with a narrow wavelength cover- 
age dTorres & R ibas 2002), which is not our case. Moreover, 
the possible systematics introduced by the use of Balmer 
lines (that can be affected by stellar winds) is compensated 
with the incorporation of He lines. In addition, the stellar 
wind would introduce a bias mainly in the systemic veloc- 
ity, which has no impact in the final distance determination. 
In any case, the derived systemic velocity (y = 164 + 5 
km s _1 ) is in excellent ag reement with the galac tic rotation 
curve of M31 (see, e.g.. iBrinks & Shangll984l where the 
studied EB is at (X, Y) (23f37Y, 8'AE)) and very similar 
to the s ystemic velocity of the neighboring EB studied in 
Paper I (y — 173 ± 4 km s _1 ). Therefore, the effect of any 
systematics on the derived radial velocities can be consid- 
ered negligible. 

- Stellar atmosphere models used to determine the surface 
flux. This is probably the major source of systematic errors. 
In order to transform the derived temperatures and gravities 
to surface fluxes, the synthetic spectra are convolved with 
filter responses and converted to abso lute magnitudes and 
color s using a series of calibrations (see lLeieune & Schaererl 
l200ll for a description on the transformations applied). The 
accuracy of these transformations is usually defin ed within 
a few hundredths of magnitude (see Bessell ^it al.ll 19981 for 
various determinations of the b olometric cor rections for the 
Sun and the colors of Vega, and iFlowerll 19961 for the disper- 
sion on the bolometric correction for hot stars). Therefore, 
we expect such systematic error to be no larger than a few 
per cent in flux, having an effect below 0.06 mag in the dis- 
tance modulus. 

- Line-of-sight absorption. The determination of the line-of- 
sight absorption benefits from the weak temperature depen- 
dence of (B - V) above ~ 30000 K to obtain an E(B - V) 
value. Therefore, the main source of uncertainty at optical 
wavelengths is the determination of the total-to-selective ex- 
tinction ratio CRv)- The value of has already been con- 
sidered to have an uncertainty of 10% and larger varia- 
tio ns are unlikely, as seen fro m previous statistical analy- 
sis (Fitzpatrick & Massa 2007). The large uncertainty in the 
value of %v increases the line-of-sight absorption error. In 
fact, uncertainty in the extinction represents ~50% of the to- 
tal error in the distance determination. Therefore, any further 
improvements that reduce the uncertainty in the line-of-sight 
absorption could potentially improve the distance determi- 
nations presented here. 

The precision and accuracy on the derived distance to M 3 1 
enables a detailed comparison with previous results. The derived 
distance modulus is in complete agreement with the mean value 
obtained f rom the combination of non-direct distance determi- 
nations in Paper II, In fact, all the 18 distances reported agree 
with the EB value within 2cr and prove that the different dis- 
tance indicators are all rather well calibrated within their error 
bars. 
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Fig. 1. Observations for M31V J004436 10+4 129 194 and corre- 
sponding W&D| fits. Top: Light curve fits and corresponding 
residuals. Bottom: RV curve fits with RVs for the primary (cir- 
cles) and secondary (squares) components with corresponding 
residuals. The phases of rejected RV observations are also indi- 
cated (crosses). 



During the course of our project, Cepheids have also been 
used to determine a distance to M 31 ( Vilar dell et al.l l2007). The 
derived distance modulus of (m - M)q = 24.32 + 0.12 mag 
represents an additional distance determination to M 3 1 that is 
fully compatible with the EB value. The consistency of this re- 
sult is specially valuable, because an EB distance to LMC has 
been used as reference. Therefore, the Cepheid distance modu- 
lus proves that the EB distances to LMC and to M 3 1 are fully 
compatible, tightening the extragalactic distance scale. In fact, 
Cepheids have proved to be extremely useful and robust to de- 
termine distances to LMC using, as reference, direct distance 
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Fig. 3. Comparison of stellar evolutionary models with derived 
physical properties of M31V J004436 10+4 129 194. Gray solid 
lines correspond to the ZAMS and the evolutionary tracks 
(ZAMS masses are indicated). The best fitting isochrone of 4.2 
Myr is also shown (gray dashed line). Gray dotted lines denote 
the uncertainties in the derived masses. The skewed rectangular 
boxes correspond to lcr error loci. 

determinations to other galaxies. Ma cri et all d2006l) used the 
Cephe ids in NGC4258, wh ich has a direct maser-based dis- 
tance dHerrnstein et al.| [T999), to determine a distance modulus 
to LMC that is almost identical to the EB value. Therefore, the 
derived distances using EBs are not only fully compatible among 
themselves, but also to other direct distance determination meth- 
ods. 

An exception to the general agreeme nt is the distance de - 
termination to M 33 with one EB system dBonanos et al. 2006). 
The derived distance modulu s is 0.3 mag larger th at the Cepheid 
distance determination of iFreedman et alj d2001l) . implying a 
distance to LMC that is incompatible with the remaining dis- 
tance determinations. Even a maser-based dista nce to M 33 was 
obtain ed favoring the short distance to M 33 dBrunthaler et al.l 
2005), the associated uncertainty is too large for a firm estab- 
lishment of the distance. Undoubtedly, additional distance deter- 
minations to M 33 with EBs will be crucial to solve the discrep- 
ancy. 

EBs have proved to be excellent distance markers and, ar- 
guably, the best method to provide direct distances for all the 
galaxies in the Local Group. In addition, the large number of 
distance indicators observed in M 3 1 strengthens the importance 
of this galaxy as the main calibrator for extragalactic distance 
determinations. In fact, al l the distance indicators used to deter- 
mine a distance to LMC dGibsonl l2000). can be used, in princi- 
ple, in M 3 1 . In addition, at least two other distance determina- 
tion methods can be used in M 3 1 (planetary nebulae luminos- 
ity function and Tully-Fisher relationship). Therefore, after the 
direct distance obtained with EBs, M 3 1 may be the best extra- 
galactic benchmark for distance determinations. 
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